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Nonlinear traveling waves that are precursors to laminar-turbulent transition and cap- 
ture the main structures of the turbulent buffer layer have recently been found to exist 
in all the canonical parallel flow geometries. We study the effect of polymer additives on 
these "exact coherent states" (ECS), in the plane Poiseuille geometry. Many key aspects 
of the turbulent drag reduction phenomenon are found, including: delay in transition to 
turbulence; drag reduction onset threshold; diameter and concentration effects. Further- 
more, examination of the ECS existence region leads to a distinct prediction, consistent 
with experiments, regarding the nature of the maximum drag reduction regime. Specifi- 
cally, at sufficiently high wall shear rates, viscoelasticity is found to completely suppress 
the normal (i.e. streamwise- vortex-dominated) dynamics of the near wall region, indi- 
cating that the maximum drag reduction regime is dominated by a distinct class of flow 
structures. 



1. Introduction 



The reduction of turbulent drag by polymer additives has received much a t tention sinc e 
it was first observed experimentally in 1940s (See reviews by 



Lumlevl <|1969) 



Virk 



IS: 



McComb (1990) 



Wei Li. 



Graham 



Li Xi and Michael D. Graham 



12004)). For a given flow rate, small polymer concentrations, 
on the order of ten parts per million by weight, can reduce the pressure drop in pipe or 
channel flow, for example, by 50% or greater. After six decades of research, the subject 
remains an active area of research, in part because of applications but also because it 
lies at the intersection of two complex and important fields, turbulence and polymer 
dynamics. A better understanding of this phenomenon may in turn yield insights into 
both the dynamics of drag-reducing fluids and of turbulent flows. The goal of the present 
work is to address turbulent drag reduction in the context of the dominant structures in 
the turbulent buffer layer, an approach which turns out to touch on many key aspects of 
the drag reduction phenomenon. 

We focus here on pressure-driven channel flow with average wall shear stress t w , of 
a fluid with dynamic viscosity rj s (in the absence of polymer), density p and kinematic 
viscosity v = T) s /p- The average streamwise velocity t/ avg and half-channel height I define 



outer scales for the flow. Inner scales are the friction velocity u T — y/ r w /p and the 
near- wall length scale l w = vju r . As usual, quantities expressed in terms of these so- 
called "wall units" are denoted with a superscript + . The friction Reynolds number 
Re T = u T ljv is simply the half channel height expressed in wall units. The Weissenberg 
number is denoted Wi = Xjw — \u^/v, where A is polymer relaxation time and j w is 
the average wall shear rate. Experimental results for a given fluid and flow geometry lie 
on curves of constant elasticity parameter El = 2Xv/l 2 . 

In channel or pipe flow, drag reduction results are often represented on a Prandtl-von 
Karman plot, U avg /u T vs. logRe T (~ log Wi — \ogEl), shown schematically in Fig. ^i. 
Point A c orresponds to trans ition to turbulence, which in Newtonian flow occurs at 



Re T ~ 45 ( Carlson et 



1982). One typical experimental path for a given polymer solu- 



tion and channel size is shown by the curve labeled "expl" . Along this path, once Re r 
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Figure 1. (a) Schematic Prandtl-von Karman plot. The dashed lines represent the experi- 
mental paths by which specific polymer systems of different molecular weights, concentrations, 
polymer-solvent pairs, etc., approach the maximum drag reduction asymptote, (b) Schematic 
of polymer induced turbulent drag reduction based on existence regions for nonlinear coherent 
states. 

exceeds a critical value (point B), the slope of the data increases from the Newtonian 
value, indicating onset of drag reduction. As Re T increases, data eventually approaches 
a new curve at point C . This curve, the so-called Maximum Drag Reduction (MDR) 
asymptote, is insensitive to polymer concentration, molecular weight or polymer species 
- all results collapse onto it at large Re T ; it is a universal feature of drag reduction by 
polymers. For large channels or low polymer concentrations, the value of Re T at the onset 
of drag reduction is independent of polymer concentration and corresponds to a critical 
Weissenberg number. For small channels or large concentr ations, how ever, diameter and 
concentration effects have been observed experimentally : specifically, there 

exists a critical pipe diameter below which, or a critical polymer concentration above 
which the flow behavior directly transits from laminar flow to the maximum drag reduc- 
tion curve as Re T increases. An experimental path showing this effect is labeled "exp2" ; 
transition from laminar flow to MDR occurs at point D. 
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Studies of drag-reducing fluids indicate that near the onset of drag reduction, t he ef- 



fects of the polymer are co nfined primarily to the buffer layer region of the flow (Virk 



1975; 



Donohue et al 



1972). Experimental observations and DNS studies show that the 



dominan t structures in the buffer layer are pairs of counter-rotating, streamwise-aligned 



vortices IjRobinsonl 1991 



Jeong et al 



1997). These vortices pull slower moving fluids away 



from the wall, forming low-speed, streamwise velocity streaks. In d rag-redu cing flows, 



these structures are modified by polymers: the buff er region thickens ( V: 



herent structures in this region shift to larger scales IjDonohue et al. 



1997; 


den Toonder et al. 


1997 


), and the burs 


ting 


Recent experimental results ( 


Warholic et al. 


1999 



1972; 



irklll975D. the co- 



Sureshku 



nar et 



I972L 



2001) reveal that in the maximum 



drag reduction region the ejections from the wall are eliminated and the near-wall vor- 
tices that sustain turbulence in a Newtonian fluid are completely destroyed. Low-speed 
streamwise velocity streaks are essentially absent. These observations suggest that the 
coherent structures in buffer layer region are crucial in addressing rheological drag re- 
duction in wall-bounded turbulent flows. 

A recent advance in the understanding of these important near-wall structures has 
come with the recognition that, in all the canonical parallel geometries (plane Cou- 
ette, plane Poiseuille, pipe) the Navier-Stokes equations supp ort nonlinear traveling wave 



state s , the family of so-called "exact coherent states" or ECS IINagata 



1986; 



Clever & Busse 



mi; Walcffe 1998112001112003: Faisst & Eckhardt 2QO.-4IWedin Kerswelll2004l). Jimenez 



and coworkers l|Jimenez fe Pinelli 



19991 



Jimenez fc Simens 



2001) have found related 



states in spatially filtered direct numerical simulations (DNS), showing the autonomous 
nature of the near- wall behavior. The flow structure of these states is a mean shear and a 
pair of staggered streamwise-aligned counter-rotating vortices, as is found in the turbu- 
lent buffer layer. In the plane Poiseuille geometry, ECS come into existence at i?e r of 44.2 



(Walcffe 
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2003), very close to the experimentally observed Re T of ~ 45 for the transition to 



turbulencel 



Carlson et al. 



1982). The spanwise wavelength L+ = 105.5 of the ECS at on- 



set closely matches the streak spacing of ~ 100 wall-units widely observed in experiments 



over a large range of Reynolds numbers (Robinson 



1991). Direct numerical simulations 



of turbulence in "minimal channel flow", i.e., flow in the smallest computational domain 
that reproduces the velocity field statistics of near-wall turbulence, give a range for the 
streamwise length L+ of 250 — 350, compared to Lt = 273.7 for the ECS, a nd a spanwise 



1991). It should be 



length that is again approximately 100 wall units IjJimenez fe M oin 
pointed out that this minimum channel contains a single wavelength of a wavy streak 
and a pair of quasi-streamwise vortices, which is the same structure seen in the ECS. 



A con ditional sampling study of coherent structures in a larger scale DNS ( Je ong et 



19971) indicates that the dominant structures near the wall in turbulent channel flow 
are counter-rotating, streamwise-aligned vortices with a streamwise length ~ 250, a 
spanwise length _L+ ~ 100 and a wall-normal size of y + ~ 50, which agrees with the 
scales of the ECS at onset. The ECS also capture the location of the pe ak, at y + ~ 12, in 



the production of turbulent kinetic energy for wall-bounded turbulence ( Kim et al 



1987 



Li et oi.ll2005D . In short, the ECS are precursors to turbulence and their structure and 
length scales closely match experimentally observed near-wall behavior. 

Because the first effects of polymer arise in the buffer region, whose structure the ECS 
evidently capture, these flows provide a natural starting point for understanding drag 
reduction. In prior work, we have studied the initia l effects of viscoelasticity on ECS in 



the plane Couette and plane Poiseui lle geometries IjStone et al. 



2003; 



Stone et al 



2004; 



Li et al 



2002; 



Stone fc Craham 



2005). The primary effect was found to be the weakening 
of the streamwise vortices, as well as changes in the statistics of the velocity fluctuations 
that are consistent with experimental observations at low levels of drag reduction. The 
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present work takes a broader view, examining the region of parameter space (Re, Wi) in 
which ECS exist and its connection to experimental observations. 



2. Formulation 

We consider pressure-driven flow with no-slip boundary conditions; v x , v y , and v z 
are streamwise, wall-normal, and spanwise components of the velocity, v, respectively. 
Reflection symmetry is imposed at the channel centerline. The laminar center line velocity, 
U, and the half-channel height, I, are used to scale velocity and position, respectively. 
The average wall shear rate "f w is given by 2U/1. Time, t, is scaled with l/U , and pressure, 
p, with pU 2 . The stress due to the polymer, t p , is nondimensionalized with the polymer 
elastic modulus, G = r/p/X, where r\ v is the polymer contribution to the zero-shear rate 
viscosity. The momentum balance and the equation of continuity are 

V-v = 0. (2.2) 

Here (3 = f] s /(r] s + ij p ) is the fraction of the total zero-shear viscosity that is due to the 
solvent, Re = pV , 1 and Re T = \p2Re. 

The polymer s tress is computed with the widely-used FENE-P constitutive model 



'Bird et al 



1987) 



a Wi (Da „ „ T \ bS 

— — + — I aVv-Vv -a =- — -, (2.3) 



1_ tra ■ 2 \Dt ' ' ' ' J 6 + 2' 

where a is a non-dimensional conformation tensor and b is proportional to the maximum 
extension of the dumbbell — tra cannot exceed b. The polymer contribution to the stress 
is given by: 

The extensibility parameter Ex — 2b ^~^ measures the relative magnitude of the poly- 
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mer and solvent contributions to the steady state extensional stress in uniaxial extension 

at high extension rate. We consider the situation 1 — [3 <C 1, in which case shear-thinning 

is negligible, as the polymer only contributes a very small amount to the total shear 

viscosity of the solution. In this situation, significant effects of the polymer on the flow 

are expected only when Ex 3> I . Finally, recall that experimental results for a given fluid 

and flow geometry lie on curves of constant elasticity parameter El = 2\v/l 2 = Wi/ Re. 

The conservation and constitutive equations are solved through a Picard iteration in 

a traveling refer ence frame - th e wave speed is part of the solution. A Newtonian ECS, 



as computed in IWaleffd l|1998j) . is first used to calculate the polymer stress tensor, t p , 
by inserting the velocity field in the evolution equation for a and integrating for a short 
length of time, usually one time unit (l/U). For this t p , a steady state of the momentum 
and continuity equations is found by Newton iteration. The resulting velocity field, v, is 
used to compute the new t p , and the process is repeated until the velocity and polymer 
field converge to a steady state. 

The momentum and continuity equations are discretized using a Fourier-Chebyshev 



FelJ 


1222; 


Min et al. 


200D 



with a third-order, compact upwind difference scheme 
the x and z directions and Chebyshev collocation in the y direction. In this as in most 
previous computational studies of polymers in turbulent flows, we have found it neces- 
sary to add an artificial stress diffusion term g^ V^, to the right-hand side of Eq. 
I|2.3|) to achieve numerical stability. The Schmidt number, Sc, which is the ratio of the 
momentum diffusivity to stress diffusivity, is set to value of 1.0. This value of Sc, though 
artificially s mall, is greater or of the same order of magnitude as that used in many 



DNS studies Sureshkumar et al 



J2002I 



1997 



Ptasinski et al 



2003; 



Sureshku mar fc Berisi l995: 



Sihilla fc Baro 



In the range of Sc where solutions can be obtained, the bifurcation 
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Figure 2. (a) Bifurcation diagram for Newtonian and viscoelastic ECS. (b) Existence bound- 
aries and drag reduction regimes for viscoelastic ECS. For all results, L x — 27r/1.0148 and 
L z = 2tt/2.633. 

diagrams shown below are insensitive to its value. The stress diffusion term is integrated 
implicitly by the Crank-Nicholson method with the other terms of the equation inte- 
grated using the Adams-Bashforth method. This equation is solved on a finer mesh than 
the momentum, continuity pair, typically 48x49x48. 



3. Results and discussion 

We study the Newtonian and viscoelastic ECS at fixed streamwise and spanwise 
lengths: L x = 2tt/1.0148 and L z = 2tt/2.633 (i.e., L+ = 273.7 and L+ = 105.5 at 
Re T = 44.2). This wavelength pair is where ECS first come into existence in the New- 
tonian case. The trivial base state in this geometry (laminar Poiseuille flow) exists at 
all Re. At Re m 977 (Re T — 44.2) for the Newtonian flow, two new solutions appear via 
a saddle-node bifurcation as shown in Fig. [2Ji. These are the ECS. These solutions are 
plotted using the maximum in the root mean square wall-normal velocity fluctuations for 
the solution, v'y . (Hereafter, an overbar indicates that the variable is averaged over the 
streamwise and spanwise directions.) The solutions with higher maximum wall-normal 
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velocity at a given Re are called "high drag" solutions due to their lower mean velocity 

at the centerline of the channel compared to the "low drag" solutions. All results in this 

paper are for the "high drag" states. Although both solutions are unstable, their status as 

precursors to transition and their structural similarity to buffer layer turbulence suggest 

that they are saddle points that underlie in part the strange attractor of turbulent flow. 

Fig. 01 indicates that the addition of polymer changes the Reynolds number i?e m i n at 

which the ECS come into existence. Curves of ECS existence boundaries i?e m i n vs. Wi 

are given for two parameter sets by the thick solid curves on Fig.JSJs. These separate the 

region where the ECS can exist (above the curves) from the region where no ECS exist, 

for the given value of Ex. While at low Wi, there is a slight decrease in i?e m ; n from the 

Newtonian value, onceW^ exceeds 45, i?e m i n for Ex = 100 is more than doubled. This 

dramatic increase in i?e m i n after onset is consistent with the experimental observation 

that the transition to turbulence in a polymer solution is delayed to higher Re than in 



the Newtonian case IjCiles fc Pettit 



1967 



White fc McEligot 



1970; 



Escudier et alWw^ 



We will refer to the Wi above which i?e m i n for the viscoelastic ECS is greater than for 
the Newtonian ECS as the onset Weissenberg number M'Zonset for drag reduction. Fig. 
|2h shows that Wionsct ~ 15, a value which is insensitive to polymer extensibility (Ex) 
or concentration (j3). Furthermore, in simulations at constant Re, it is found that the 
value where the centerline mean velocity Z7 max first exceeds the Newtonian value - the 
nearly vertical set of points in Fig. 03 - is also located at Wi « 15 in the parameter 
regime that has been examined here. This onset value is high by ab out a factor of two 



compared to values pr edicted by two recent viscoelastic DNS studies IjHousiadas fc Beri? 



2003 



Min et al, 



2003), but in those studies El was significantly smaller, and the onset 
Reynolds number correspondingly larger, than the values considered here. 

Figure |2t shows mean velocity profiles at six different sets of parameter values, each 
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Figure 3. Mean streamwise velocity for Newtonian and viscoelastic ECS on the ECS 

existence boundary. 

corresponding to a point on the existence boundary for the ECS (i.e. a bifurcation point). 
Remarkably, they all fall on the same curve, when plotted in outer units. Therefore, at 
least for the values of Re and Wi that are currently accessible in our simulations, we 
observe that mean velocity profiles at onset of the ECS have a universal form. 

We now turn to the study of the evolution of the ECS along some experimental paths, 
lines of constant El. Two such paths, the thin solid lines with hollow symbols, are shown 
in Fig. 03. Consider first the case El = 0.010; as Re and Wi increase, the path inter- 
sects the ECS existence boundary at point A and the drag reduction onset threshold 
curve at point B, where the transition to turbulence and the onset of drag reduction 
occur, respectively. Turning to the case El = 0.019, mean velocity profiles expressed in 
wall units are shown for various values of Re in Fig. |2t>- For this parameter set, drag 
reduction is observed immediately upon onset of the ECS. Along with drag reduction, 
enhanced streamwise velocity fluctuations and the reduced wall-normal and spanwise ve- 
locity fluctuations are found, consistent with experimental observations and DNS results 

fri in 

at low to moderate degrees of drag reduction (Virk 1975; Surcshku mar et oiJH997(l . The 



effect of viscoelasticity can also be observed in the reduced Reynolds shear stress and 
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Figure 4. (a) Streamwise velocity for an exact coherent state at Re = 2200(7?e T = 66.3), 
Wi = 41.8, Ex = 100, (3 = 0.97. Range: (black) - 0.58 (white), (b) Trace of the polymer stress 
for the same state. Range: (black) -1800 (white). 



ultima tely can be traced to th e 



ticity IjStone et al 



2002, 



2004; 



suppression o f the streamwise vortices by the viscoelas- 



Li et al 



2005j) . Figure ^ shows fields of v x and trx p on 



the El = 0.019 path at Wi = 41.8, Re = 2200 (the point just left of the label "C" on 
Figure Et). The region of high polymer stress clearly "wraps around" the streamwise 
vortices, and the corresponding polymer force (~ V • r p ) is in direct opposition to the 
vortex motions. 

Continuing upward in Re and Wi at El — 0.019, the path re-intersects the ECS ex- 
istence boundary, at point C in Fig. |5Jj. (We suspect that this will also happen in the 
El = 0.010 case, but at higher Re and Wi than are accessible with our current com- 
putational approach.) Above this point the flow can no longer sustain these ECS; vis- 
coelasticity completely suppresses the near-wall vortical structures. This result is con- 
sistent with experimental observations that, in the MDR regime, the eruptions of low- 
momentum fluid from the wall arc eliminated and the near-wall streamwise vortices are 



completely destroyed IjWarholic et al. 



1999 



20011) . Experimental results also show that 



in the MDR regime, the Reynolds shear stress is much smaller than the Newtonian 



12 



value (Warholic et al. 



We 



1999 



-fer 



| AZ uud MiUimil D 



2001 



Graham 



Ptasinski et aZ. 2001), and stre amwise velocity fluctu 



1999) 



ations decrease to levels close to or below the Newtonian value \ Warholic et 
All these observations suggest that the turbulent production and dissipation take place 
by a different mechanism in the MDR regime than at lower degrees of drag reduction. 
Although our study does not reveal this mechanism directly, it does suggest that the 
disappearance of ECS is related to the MDR regime. This result encourages us to take a 
broader view, examining the region of parameter space (Re, Wi) in which ECS exist and 
its connection to experimental observations. 

Fig.^3 is a schematic based on the results shown in Fig. ^jp. Line 1 represents the ECS 
existence boundary at constant Ex. Line 2 represents the drag reduction onset threshold, 
which separates the ECS existence region into "turbulence without drag reduction" and 
"turbulence with low or moderate drag reduction" regions. Line "expl" represents an 
experimental path at constant El, which passes through the ECS existence region. In this 
case, as Re (and Wi) increases, this path intersects with the ECS existence boundary at 
point A and drag reduction onset threshold at point B, where the transition to turbulence 
and the drag reduction onset occur, respectively. Note the correspondence with points A 
and B on the schematic Prandtl-von Karman plot, Fig.^,, as well as on Fig. 03. As Re and 
Wi continue to increase along this path, the system will eventually exit the ECS existence 
region at point C, where the flow can no longer sustain these ECS. Experimental results 
show that in the MDR regime, near-wall streamwise vortical structures are essentially 
absent. Our results together with this experimental observation suggest that the loss 
of ECS may be somehow related to the approach of the MDR regime, in which other 
types of coherent traveling wave states (temporally intermittent structures, hairpins, 
Tollmien-Schlichting waves, intrinsically elastic structures,...) may be unmasked and 
become dominant. This possibility is represented by line 3 in Fig. \\]p, a hypothetical 
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existence boundary for a distinct class of flow structures that exists at high Wi. In 

this scenario, the crossing of path expl across point C represents the transition to the 

maximum drag reduction regime. This scenario, incorporating transition to turbulence, 

onset of drag reduction and approach of the MDR regime is consistent with the behavior 

on experimental path expl shown in Fig. ^i. 

Now consider the experimental path "exp2" on Fig. [flp. This path corresponds to a 
value of El that does not intersect with ECS existence region at all. For the conditions 
/3 — 0.97, Ex — 100 shown on Fig. Et: this situation arises if El > 0.024. The scenario on 
Fig-d 3 would predict in this case that, with the increase of Re and Wi, the flow behavior 
directly transits from laminar to MDR at point D. As El is inversely proportional to I 2 
(or R 2 in pipe flow) this prediction is consistent with experiments in small diameter pipes 
- the "diameter effect" as exemplified by experimental path exp2 in Fig. ^i. 

The "concentration effect" can also be captured by this scenario, as we now describe. The 
quantity S — 1 — (3 is proportional to polymer concentration in dilute solution. Using S, 
the parameters Ex and El can be written as Ex = ^bS/ (1 — S) and El = 2\n s / pl 2 (l — S). 
Thus while El is virtually unchanged by a change in S, Ex is proportional to it. An 
increase in Ex compresses the ECS existence boundary leftward and eventually a given 
experimental path can no longer intersect the ECS existence region, resulting again in 
flow behavior that directly transits from laminar to MDR. 

Finally, we observe that the existence boundaries can be interpreted in terms of length 
scales. Recall that the half-height of the channel, expressed in wall units, is simply Re^- = 
y2Re. Thus the existence boundary corresponds to the minimum half-channel height in 
which an ECS can exist, as a function of Wi. Points where a line of constant El intersects 
the existence boundary are points where the channel height and the minimum height for 
the existence of an ECS coincide. 
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4. Conclusions 

Many observations of drag reduction in dilute polymer solutions are m irrored by the 



effect of viscoelasticity on the channel flow ECS discovered by Waleffe (jWaleffe 



2001, 



20031) . The transition behaviors from laminar to turbulent flow, from no drag reduction 
to drag reduction, and from moderate drag reduction to MDR can be connected to the 
birth, evolution and death of these ECS, respectively. Our results and the scenario that 
we infer from them yield explicit predictions, testable by DNS, with regard to all these 
phenomena. 
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